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We comment on the recent paper (hep-lat/0102003) by Horvath, Isgur, McCune, and Thacker 



which concludes that the local chiral structure of fermionic eigenmodes is not consistent with in- 
stanton dominance. Our calculations, done with an overlap action, suggest the opposite conclusion. 
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]^ • I. INTRODUCTION 

o: 

' In a recent paper, Ref. Horvath, Isgur, McCune, and Thacker present evidence from lattice simulations claiming 
to show that the behavior of fermionic eigenmodes in the QCD vacuum is inconsistent with an instanton picture. They 
refer to a recent paper of ours, Ref. 1^, and remark "Although their qualitative conclusions regarding the instantons 
differ from ours, there does not seem to be any direct conflict between their lattice results and ours." We disagree 
with that statement. To expose that disagreement, we reanalyzed a subset of our data along the lines of Ref. 

Why is it important to resolve this disagreement? There is an apparent inconsistency, first described by Witten in 
. 1979 1^, between large number of colors (large- iVc) dynamics and instanton-based phenomenology. We do not wish 
^ ' to summarize arguments favoring one or the other approach, which are well-described in a variety of publications 
CnJ 0- In principle, lattice simulations might distinguish between the qualitative features of the two phenomenologies. 
However, the lattice calculations of Refs. [Q and |^ produce contradictory results. 
, The two simulations differ in several ways. Both are done in quenched approximation, using the Wilson gauge 

■ action. Ref. ||l| works at a lattice coupling (3 — 5.7 and Ref. |^ works aX P — 5.9. The lattice spacings a are at 
nominal values of 0.17 and 0.11 fm, respectively, from the Sommer parameter, using the interpolating formula of Ref. 

. 1^. The fermion action of Ref. 0| is the Wilson fermion action. This action has rather poor chiral properties, as 
' evidenced by its large additive quark mass renormalization and the presence of exceptional configurations. Formally, 
. it has 0{a) scaling violations. Its real eigenmodes are not eigenfunctions of 75, and the expectation value of 75 of a 
' I ' real eigenmode is a function of the eigenvalue. 

Oh! The action of Ref. Q is an overlap action ||^ , which is constructed by using an approximate overlap action as its 
^ ■ starting point . This action has exact chiral symmetry. There is no additive mass renormalization of the quark mass 
and there are no exceptional configurations. Zero-eigenvalue eigenmodes are chiral eigenstates. Nonzero eigenvalue 
eigenmodes have zero expectation value of 75. The action has only 0{a^) discretization artifacts. 

Our action uses "fat links." Fat link actions replace the usual one-link gauge connection with a combination of 
several gauge paths. This combination might or might not be projected to SU(3), but in any case it is gauge invariant 

■ and local. Locality assures that the fat link action is in the same universality class as the thin link one. Fat links 
improve the chiral behavior of nonchiral actions Q, and their use in our implementation of the overlap action is 
simply to make the calculation of the overlap operator more efficient (a factor of roughly 20 for the action of Ref. 
0, compared to the use of the thin link Wilson action in the overlap.) The particular choice of fattening is not too 
important. In Ref. we used an APE-blocked link Q as the gauge connection. In Ref. ||l^ it is shown that the 
perturbative effect of N levels of APE-smearing with a smearing parameter a is to multiply the vertices by a form 
factor (1 — with q the gluon momentum and — 4/a sin^ In coordinate space this corresponds 
to a Gaussian smearing with a spread (x^) = aN/3a^ of the quark gluon vertex. Some readers might be concerned 
that several levels of APE smearing might adversely affect our results. To address that point, we describe results 
from another fattening, one chosen to lie rigorously within a hypercube | pd| . We observe no qualitative change in our 
results. 



II. ANALYSIS OF FERMION EIGENMODES 

The statement of the problem has been given by Ref. |jl| : "If an instanton-dominated picture of low eigenmodes is at 
all valid, we would expect the peaks of the wave function to closely resemble instanton zero modes. Thus, if instantons 
dominate, a local peak in the wave function for a low-lying eigenmode (a fermion lump) should be dominantly a lump 



1 



in ip\ipL or ipl^ipR, but not both." (Here V'lV'i = i^H^ ~ 75)V' and ipl^ipR = i^Hl + 75)^'-) "On the other hand, 
fermion lumps without a definite chirahty would be an indication of non-self-dual fluctuations, as expected from a 
confinement-related mechanism." 

In Ref. |§|we used an overlap action with = 10 levels of APE smearing with a — 0.45 APE parameter. We 
measured the autocorrelation function of the local chirality density lu{x) = (^{x)\"f5\^{x)) and the correlator of uj{x) 
with a pure gauge observable sensitive to topological charge Q{x). We saw that both kinds of correlators were strongly 
peaked at small operator separation indicating the the chiral density has spatially localized lumps and those lumps 
correlate with the lumps of the topological charge density. However, there is a logical possibility that a lump in 
Lu{x) = V'lV'l — ipRipR does not indicate a lump in either i/'li/'l or tpRtpR- To check that possibility, we follow Ref. § 
and define a "chiral order parameter" X{x) 



tan(-(l 



-X{x)) 



\'^Fiix)\ Vfl(a;)V'i?.(a; 



(1) 



and plot histograms of the density of the variable X, p{X), for regions of the lattice where ■ilj{xyip{x) is large, a few 
per cent of the points of the lattice. We distinguish two cases: 

(a) Zero modes are chiral. The analogs of Figs. 1 and 4 of Ref. 0| for our action areHust delta functions p{X) ~ 
5{X ±1). The chiral symmetry breaking inherent in the Wilson fermion action of Ref. |l| broadens this distribution. 

(b) Nonzero eigenvalue eigenmodes are nonchiral, but we find that the distribution p{X) is strongly peaked near 
X = ±1. To show this, we take a set of six 12'' lattices from our data set and plot histograms of p{X). Fig. [l] shows 
the result for the lowest two nonchiral eigenmodes for the overlap action, with a cut keeping the top 2.5 per cent of 

The imaginary part of the eigenvalue of the Dirac operator for these modes is in the range 0.03/a to 0.09/a. 
Keeping more lattice points in the histogram will fill in the dip. That simply implies that as we keep more and more 
lattice points we pick up more of the vacuum fluctuations around the chiral lumps and those eventually overwhelm 
the chiral structure. This happens very slowly with the two smallest non-chiral modes. Even with 30 per cent of the 
lattice points the two peak structure is clearly visible. 
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FIG. 1. Histogram of p{X) of the lowest two nonchiral eigenmodes for the overlap action, with a cut keeping the top 2.5 
per cent of i/)^^/). 

Figs. H and || shows the result for the largest modes we have recorded, the ninth and tenth eigenvalue modes of 
the squared Dirac operator. The imaginary part of the eigenvalue of the Dirac operator for these modes is around 
0.20/a-0.26/a. Figure ||, where about the top 0.3 per cent of ■i/'^V' is kept in the histogram, shows strong peaks around 
X = ±1. The peaks disappear as we add more points to the distribution. With 2.5 per cent of the lattice points kept 
the histogram is almost flat as Figure |^ shows. 

Lower eigenvalue eigenmodes show stronger peaking than higher eigenvalue modes, but all the modes we examined 
show peaking. We remind the reader that we showed in Ref. pj that hadronic correlators for light quark masses are 
well approximated by quark propagators containing only a few eigenmodes. 
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FIG. 2. Histogram of p{X) of the highest two eigenmodes (out of ten) for the overlap action, with a cut keeping the top 
0.3 per cent of tp'^^p. 
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FIG. 3. Histogram of p{X) of the highest two eigenmodes (out of ten) for the overlap action, with a cut keeping the top 
2.5 per cent of 



Next we consider results obtained with an overlap action that uses more local fat links. The fat links of the 
hypercubic blocking plf mix links only within a hypercube but achieve almost the same level of smoothness as the 
APE smearing considered previously. 

The analog of Fig. |l| from eigenmodes of the overlap action with hypercubic fat link is shown in Fig. |[ We used 
the same 12* configurations and kept the same fraction of lattice points as in Fig. |l|. The shape of the curve is 
qualitatively similar to that from APE-blocked links. Peaks of fermion density are again chiral. 
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FIG. 4. Histogram of p{X) of the lowest two nonchiral eigenmodes (out of ten) for the overlap action, with a hypercubic 
fat link, with a cut keeping the top 2.5 per cent of tli^ip. 



III. CONCLUSIONS 



Our studies with a chiral fermion action at a lattice spacing near 0.11 fm showed that low eigenvalue fermionic 
eigenmodes have a structure strongly correlated with the locations of instantons and anti-instantons. This note 
demonstrates that in the locations where the fermionic modes are largest, the modes are also chiral. It seems to us 
that the simplest description of what we see is in terms of an instanton liquid model of the QCD vacuum. 

We do not doubt that the numerical results of Ref. [|l| arc correct, but we are concerned that the combination 
of a larger lattice spacing and the use of a nonchiral fermion action introduces lattice artifacts which compromise 
extrapolation to the continuum limit. 
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